CSPIN/0501 



Inelastic transport in molecular spin valves 
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We present a study of the effects of inelastic scattering on the transport properties of various 
nanoscale devices, namely H2 molecules sandwiched between Pt contacts, and a spin-valve made by 
an organic molecule attached to model half-metal ferromagnetic current/voltage probes. In both 
cases we use a tight-binding Su-Schrieffer-Heeger Hamiltonian and the inelastic effects are treated 
with a multi-channel method, including Pauli exclusion principle. In the case of the H2 molecule, 
we find that inelastic backscattering is responsible for the drop of the differential conductance at 
biases larger than the excitation energy of the lower of the molecular phonon modes. In the case of 
the spin-valve, we investigate the different spin-currents and the magnetoresistance as a function of 
the position of the Fermi level with respect to the spin-polarized band edges. In general inelastic 
scattering reduces the spin-polarization of the current and consequently the magnetoresistance. 

PACS numbers: 



I. INTRODUCTION 

Recent years have witnessed an increasing interest 
for spintronics devices and for their possible appli- 
cation as novel elements in microelectronics industry 
S Read-heads based on the giant magnetoresis- 
tance (GMR) effect 0, |5j have already revolutionized the 
hard-drives market and other magnetic devices such as 
magnetic random access memory (MRAM) are expected 
to have a similar or even bigger impact. These devices 
exploit both the magnetic and the electronic properties 
of conducting materials. 

In the typical spin-valve geometry a non-magnetic 
spacer is sandwiched between two magnetic contacts. 
These provide a source of spin-polarized electrons and the 
device practically behaves as a polarizer/analyzer spin 
filter, with the magnetization vector establishing the po- 
larization axis. Typically the resistance of the device 
is higher when the magnetization vectors of the contacts 
are aligned antiparallel to each other, and drops when the 
alignment is parallel. The non-magnetic spacer controls 
the total resistance of the device and the overall trans- 
port regime (diffusive, ballistic, tunneling, etc.), while 
the quality and nature of the interfaces can further affect 
the spin polarization. 

Many suggestions for possible non-magnetic spacers 
have been recently proposed Q. Among the most in- 
triguing there is the possibility of using organic molecules 
0- These offer low fabrication costs and high manufac- 
turability typical of molecular electronics H,E1 anci the 
advantage of long spin-relaxation times. Several proto- 
types for molecular spin-valves have been already fabri- 
cated and GMR in excess of 10% have been reported for 
spacers such as 7r-conjugated molecules 0, an( i m 
organic tunneling junctions 12] . 

Electron conduction through these systems is always 
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explained in terms of elastic transport 0, ll3L fl4| and in- 
elastic effects are often neglected. This may be a quite 
drastic simplification in particular in the case of thick 
organic spacers, for which the transmission due to di- 
rect tunneling becomes negligible. Inelastic scattering 
not only does open new transport channels but also af- 
fects the transmission of the elastic ones by introducing 
competition for the same final states due to the Pauli ex- 
clusion principle [l5|, [ljj . These effects can change dras- 
tically the conducting properties of the devices and affect 
the GMR of spin-valves. Recent experimental and theo- 
retical works have already shown the importance of the 
inelastic scatte ring in different non-magnetic molecular 
devices [II El EIM El El El El 

In this work we investigate the effects of inelastic scat- 
tering on the transport properties of model spin-valves. 
First, we briefly introduce our method for calculating 
the I-V characteristics of a typical two probe device in 
presence of electron-phonon coupling. This is based on 
a combination of scattering theory in a Green's function 
framework El E3 and a self-consistent evaluation of the 
non-equilibrium electron distribution |16| . Then we show 
a typical calculation by investigating the transport prop- 
erties of the H2 molecule sandwiched between two plat- 
inum leads. Finally we consider the case of organic spin- 
valves obtained by sandwiching a mono-atomic chain, re- 
sembling the properties of a (CH)„ molecule, between 
half-metallic leads. For this system we consider a model 
Mpi half metal (according to the classification of Coey 
and Sanvito E3)j which resembles the electronic struc- 
ture of magnetite. 



II. COMPUTATIONAL METHOD 

A. Problem set-up 

The system under investigation is schematically rep- 
resented in figure ^ It consists of two semi-infinite 
current/ voltage probes formed from mono- atomic chains 
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and sandwiching a scattering region (the molecule). 
Electron-phonon interaction is present only in the scat- 
tering region while the leads act as an ideal phonon sink. 
The scattering region is composed by one or more atoms 

|k > #-#-#WWW#-#-# lk>> 

FIG. 1: Schematic representation of a system composed of 
two semi-infinite current /voltage leads and a scattering re- 
gion, where inelastic scattering can occur. An electron with 
fc-vector k approaching the scattering region from the left is 
transmitted to the right with k- vector k' . In this ID case if 
k ^ k' the energy is not conserved. 

connected to the leads both via electron and phonon- 
assisted hopping. 

The Hamiltonian H describing such a system can be 
divided in five parts 

H = H L + Vml + H M + V M r + H R . (1) 

Hjj (Hji) describes the left-hand side (right-hand side) 
lead, Hm describes the molecule and Vml (Vmr) the in- 
teraction between the left (right) lead and the molecule. 
For the leads we consider a tight-binding model expanded 
over a s-like orthogonal basis set \i), where i labels the 
atomic position. Although there is no conceptual compli- 
cation in considering a multi-orbital tight-binding model 
here we use only one orbital per atomic site in order to 
keep the computational overheads to a minimum. This 
choice is motivated by the severe scaling of our compu- 
tational method |l6j | with the number of phonons. Thus 
7?l ^d i?R read 

-1 

Hl = e L J2 c ^ ( c I c *+i + c ' c -) . ( 2 ) 

%— — oo i— — oo 

+ 00 +OG 

i=N+l i=N+l 

In the expressions above c, and c\ are the annihilation 
and creation operators for electrons at the site i, (3 a (a= 
R, L) is the hopping parameter in the leads and e a is the 
lead on-site energy. This already include the rigid shift 
given by an external voltage V applied to the system, 

that is €l = £l + T an d e R = £ r — y> w ^ n e L/R tne 
on-site energy for V = 0. Note that we implicitly assume 
that the leads are good metals presenting local charge 
neutrality 28] . The leads extend from — oo to i = to 
the left and from N + 1 to +oo to the right, with the 
molecule comprising N sites (i = 1,N). 

In our model system the electron-phonon interaction is 
present only in the molecule and at the interface of this 
with the leads. We use for the description of such a scat- 
tering region the Su-Schrieffer-Heeger (SSH) Hamiltonian 
|29l l3Ct 1.3 1| in its second quantization form, considering 



only longitudinal phonon modes |2J|. Thus Hm reads 

N 

H M = em X! b l bi + X! *y ( b i b J + b ) bj ^ + X! ^fe a l a fe + 

i— 1 i^j k 

k i=£j 

Such a SSH-like Hamiltonian describes linear electron- 
phonon coupling. It has a purely electronic part 
written in a tight-binding form with annihilation and 
creation operators hi and b\, hopping parameter fy 
and on-site energy cm- These tight-binding parame- 
ters refer to the static system (without phonons). In 
addition there is an electron-phonon interaction part 
J2k lij{ a \+ a k)(b\bj+bjbi), which describes phonon 
absorption/emission. This is linear in the phonon annihi- 
lation and creation operators a& and at with the coupling 
strength given by 7,^. Finally it comprises the purely 

phononic energy ^ fc tvjjka\,ak with ojk the phonon fre- 
quency for the fc-mode. The frequencies Wk are calculated 
using a classical "springs and balls" model for the scat- 
tering region |32j |. i.e. by diagonalizing the dynamical 
matrix obtained with the appropriate boundary condi- 
tions (see figure [3). 




FIG. 2: Classical "springs and balls" model used to calculate 
the phonon modes. The leads act as hard walls constraining 
the oscillation to the middle atomic chain, m-j are the masses 
of the atoms in the scattering region and Kj are the spring 
constants. 

The parameter 7^ depends on both the energy of the 
/c-mode, Huik, and the vibrational properties of the atomic 
sites i and j. It can be written as 

7*=aJ^L-^UyA (5) 

where ei is the displacement vector of the i-th atom asso- 
ciated to the fc-mode, m., (m,-) is the mass of the i-th (j- 
th) atom and aty is the coefficient of the linear expansion 
of the hopping parameter over the atomic displacement 

Kli ~ 1j) = l ij + a ij (li ■ ( 6 ) 

In the previous equation Uj is the first term of the expan- 
sion and stands for the hopping parameter of the static 
system at the equilibrium distance, while qi and qj are 
the displacements from the equilibrium positions of the 
atoms i and j due to the lattice vibrations. 
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The Hamiltonian coupling the molecule to the leads is 
again of a SSH form 



Vml = 



, + ]Tr fe (4 + afe ) 



(cj&i + b\c Q ) , (7) 




with an analogous expression for Vmr- The purely elas- 
tic part is described by the hopping parameter w a while 
the inelastic part is characterized by the electron phonon 
coupling T k 



(8) 



The equation JHJ is obtained from the (JHJ) by imposing 
static boundary conditions to the leads (qi = for i be- 
longing to the leads). 



B. Inelastic scattering as an effective elastic 
problem 

We tackle the problem of transport in presence of 
electron-phonon interaction by usin g th e method pro- 
posed first by Bonca and Trugman |33l 0|. This ap- 
proach consists in mapping a general many-body trans- 
port problem onto a fictitious one-body problem, where 
only elastic scattering is present. The method, in con- 
trast to other perturbative approaches [35| . is in princi- 
ple exact, i.e. it does not introduce any approximation 
in calculating the inelastic contributions to the current. 
The drawback is cast in the dimension of the Hilbert 
space needed for defining the fictitious one-body system, 
which grows rapidly with the number of phonon states 
considered. In this section we will restrain our attention 
to mono-atomic single-channel leads, but the same pro- 
cedure can be easily extended to the multichannel case 

The main idea behind this approach is that the energy 
is globally conserved during the scattering process. In 
fact, although an electron leaving the left-hand side lead 
with energy E and emitting a phonon is absorbed by the 
right-hand side lead with a different energy E', the energy 
of the system electron plus phonon is conserved (E—E' is 
the phonon energy) . This suggests the possibility of con- 
structing an associated multi-channel scattering problem 
[3?| , using asymptotic states that include both electronic 
and phononic degrees of freedom, all having the same to- 
tal energy. In this way the scattering problem is formally 
elastic and can be solved by standard scattering theory 
[T5I l25| . However, when extracting the electron current, 
one should be able to map this fictitious problem back to 
the original many-body one. 

We then define states \k; n a> np, ...) where the elec- 
tronic part \k) describes an asymptotic Bloch state k in 
the leads, and the phononic part |n a , np, ...) is a Fock 
state associated to the various phonon modes a, (3, ... 
of the molecule. Here n a is the occupation number of 



the mode a. During a generic scattering event the to- 
tal energy of such state E = Ek + h^ a (n a w a ) is con- 
served, despite the fact that the electronic energy Ek may 
not be. Therefore the scattering between the asymptotic 
states \k;n a ,rip, ...) of this fictitious system is formally 
elastic. Note however that the Bloch states \k) asso- 
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inelastic 



E\ lk'> 



elastic 
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FIG. 3: Mapping procedure for a mono-atomic system with 
only one phonon mode, (a) During an inelastic process elec- 
trons change their energy from E to E' . In contrast in our 
fictitious system the scattering between states with different 
electronic energy (\k, 0) — > \k' , 1}) is elastic, since the total en- 
ergy (electron plus phonon) is conserved, (b) A mono-atomic 
chain sandwiching a molecule with vibrational degrees of free- 
dom is mapped onto an equivalent system of disconnected 
atomic chains coupled through a scattering region. 

ciated to the different asymptotic states \k;n a ,np, ...) 
in general have different energies, and must be calcu- 
lated from the electronic secular equation at the energy 
Ek = E — h^2 a (n a Lu a ). Importantly if the scattering is 
between states having identical phononic parts, then it is 
elastic not only for the fictitious system, but also for the 
actual many-body problem. 

Therefore the mapping |fe) — > \k;n a ,np, ...) transforms 
an inelastic problem onto a fictitious multi-channel elas- 
tic one (see figure |3J). For example within this approach 
a linear atomic chain sandwiching a vibrating molecule is 
formally equivalent to a collection of linear chains sand- 
wiching an elastic scatterer, where the coupling between 
the chains is provided by the electron-phonon interaction 
only (see figure EJ. 

One can now construct a scattering theory for the 
asymptotic states \k; n a , np, ...) by using one of the sev- 
eral methods available for elastic transport [2(|. In this 
case however there are two additional features. First a 
new expression for the current, including Pauli exclusion 
principle, should be derived. In fact different incoming 
asymptotic states corresponding to electronic states with 
different energy may compete for the same final Bloch 
state. Secondly the dimension of the Hilbert space is in 
principle infinite, since the occupation numbers entering 
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the states \k; n a , np, ...) are unbound. In practice how- 
ever only a finite number of phonon modes is usually 
sufficient to obtain converged results, that is when the 
scattering transmission coefficients do not vary with the 
number of phonon states. Most importantly the conver- 
gence can be carefully monitored by progressively enlarg- 
ing the Hilbert space. 



C. Scattering theory 

In the fictitious system, where the asymptotic states 
are \k;n a ,np, ...), the scattering is formally only elastic. 
Therefore we can calculate the quantum mechanical scat- 
tering probabilities by using standard scattering theory. 
In this work we have used the Green's function approach 
introduced by Sanvito et al. [2{|. The method consists 
in evaluating the total retarded Green's function for the 
entire system (leads plus scattering region), from which 
extracting the scattering matrix. Here we briefly sum- 
marize the method and highlight the main modifications 
needed when dealing with the fictitious system, instead 
of a purely electronic one. Furthermore we consider sim- 
ple one-dimensional leads, supporting only one scattering 
channel at a given energy. 

The starting point consists in writing the surface 
Green's function for the leads. These are in the form 
of a block diagonal matrix, which, in the case of the left- 
hand side lead, reads 



Gl(E) 



(9l{E) 

v ... . 



9l 



(E") 



J 



An analogous expression is valid for the right-hand side 
lead surface Green's function G R (E). g^(E') is the left 
surface Green's function of the purely electronic system 
(no phonons) at the energy E' . In general it is a N x N 
matrix, with N the number of degrees of freedom (or- 
bitals) describing the leads surface, but in this simple 
case it is only a c-number 



(9) 



with k the adimensional fc-vector solution of the secular 
equation for the simple chain 



E' = £l + 2/3 L cos(fc) 



(10) 



[i is a collective index, which label the phononic 
state [i = {n a ,np, ...} and the energy E' is the elec- 
tronic energy of the state \k; n a , np, ...), i.e. E' = 
E — h^ a {n a oj a ). Note that the block diagonal form of 
G*l and G R is the result of the fact that different phonon 
states \n a ,np, ..) do not mix in the leads. Note also that 
the Green's functions G° and G R are in principle infi- 
nite matrices since the number of possible phonon states 



is infinite. In practice states with very large occupation 
numbers are not accessible and a cutoff N p h over the 
number of phonon states is used. 

In contrast to what happens in the leads, different 
phononic states do interact inside the scattering re- 
gion. The matrix representing such an interaction is a 
(M • A^ph) x (M ■ N v h) hermitian matrix, where M is the 
number of electronic degrees of freedom (orbitals) de- 
scribing the molecule. However, it is important to note 
that in this simple scattering approach, where no self- 
consistent evaluation of the scattering potential is car- 
ried out, it is possible to eliminate the explicit depen- 
dence of the scattering region Hamiltonian over the inter- 
nal electronic degrees of freedom. This is done by using 
the recursive decimation method [2j|, to yield an effec- 
tive energy-dependent 2iV p h x 2 N p h Hamiltonian H c h (E), 
coupling all the phononic states in the left-hand side lead 
with those of the right-hand side lead 



H „(E)-( H£ R (E)\ 



(11) 



The blocks and H R describe the interaction respec- 
tively within the last plane of the left-hand side and first 
plane of the right-hand side lead, while the off diagonal 
blocks describe the coupling between these planes. 

Finally the total Green's function of the system leads 
plus molecule calculated at the two extremal planes of 
the leads (i=0 and i=N+l) can be obtained by simply 
solving Dyson's equation 



G{E) = [G\E)- 1 - HcsiE)]- 1 =( ^ ^ LR 



> Grl Gr 

where G°(E) is written as 

o _ ( Gl(E) \ 
G {h) - \ G° R (E) ) ■ 



(12) 
(13) 



The total Green's function G contains all the informa- 
tion about the scattering amplitudes. These can be ex- 
plicitly extracted by using the Fisher and Lee relation j3^ 
or in more general terms by proj ecting G over a general 
scattering wave- function |25l l26|. For our single-mode 
leads these simply read 



-6p V + ih(G^)^^/v~v^ 



r' = -6^ + ihiG^^y/Vf^Vv 



(14) 



(15) 



(16) 



(17) 
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where is the group velocity associated to the elec- 
tronic component of the state = \k; n a , np, ...) 



(18) 



^„ and r M „ are respectively the transmission and re- 
flection amplitude for states incoming from the left-hand 
side lead, while t' and r' are the same quantities but 
for states incoming from the right-hand side lead. 

Finally the current can be extracted from the Landauer 
formula 0] 



2e 
~h 



J ^2,[\t^u\ 2 F L {E) - F R (E)]dE 



(19) 



where Fl(E) and Fr(E) are the Fermi distributions re- 
spectively for left- and right-hand side lead. Note that 
in the case of time-reversal symmetry t^u = 
current becomes 



t' and the 



2e 

7T 



jY^\t,AE)\ 2 [F L {E)~F R {E)]dE. 



(20) 



D. Pauli principle 

For elastic transport different scattering channels do 
not compete for the same final state, since the scattering 
wave-functions extend from one lead to the other. This 
is also formally true for our fictitious system. However 
when one considers the electronic part of the fictitious 
system, the situation appears rather different. Consider 
for example the scattering process in which \k;n a ,np, ...) 
is transmitted into the same state \k; n a , np, ...), and the 
one in which it is the state \k'; m a , mp, ...) to be transmit- 
ted into \k; n a , np, ...}. In the first process the phononic 
configurations of the initial and final states are identical 
and therefore the scattering is elastic. In contrast the sec- 
ond process involves a change in phononic configuration, 
which means that energy has been exchanged between 
the electronic and the phononic sub-systems. This is an 
inelastic scattering event. 

The crucial point is that these two processes lead to the 
same final electronic state \k), the first through elastic 
and the second through inelastic scattering. These two 
scattering events do not share the same scattering wave- 
function and therefore compete for the same electronic 
states. Therefore in this case the Pauli exclusion prin- 
ciple should be taken explicitly into consideration when 
evaluating the current. 

We then use the procedure introduced by Emberly 
and Kirczenow [16j , which evaluates self-consistently the 
non-equilibrium electron distributions resulting from the 
Pauli principle. We define {f" a e i) R (E') as the electronic 
distribution of the right-hand side lead resulting from 
electrons transmitted elastically from left to right at an 
energy E'. As a matter of notation we label with "+" 



("— ") scattering processes involving electron transmis- 
sion (reflection), "R" ("L") indicates that the distribu- 
tion is for the right-hand side (left-hand side) lead, and 
"el" that the scattering process is elastic. Finally the 
indexes a label the associated states |a) = \k; n a , np, ...) 
describing the process. Similarly (f" a in ) R {E') is the elec- 



a) 



lk> I0> 
lk'>ll> 



lk> I0> 
lk'>ll> 



b) 




FIG. 4: Scattering events of the fictitious system, contribut- 
ing to the various electronic distributions, (a) elastic scat- 
tering of the state |fc;0) = |fc)|0), contributing to the dis- 
tribution (f" a e i) R (E'). (b) inelastic scattering of the state 
|fc;0) = |fc)|0) into the state \k ; 1) = contributing to 

the distribution {f+ a L) R (E'). 

tronic distributions in the right-hand side lead given by 
electrons incoming from the left-hand side lead and in- 
elastically transmitted through the molecule (Fig^J. In 
this case E' is the energy of the final electronic state and 
the scattering is between the states |a) = \k;n a ,np, ...) 
and | a') = \k' \m a ,mp, ...) . 

If one neglects Pauli's principle these quantities are 
simply related to the Fermi distributions of the lead con- 
taining the incoming electrons 

{f+%) R {E') = F L (E')T aa (E',E') , (21) 



{n%) R {E') = F L {E)T aa ' (E,E'), (22) 

where we have introduced the transmission coefficients 
T af3 (E, E') = \t a p\ 2 , which depend on both the initial E 
and final E' electronic energy. All the other distributions, 
for electrons approaching from the right-hand side lead, 
can be written in a completely analogous way 

We now need to establish a procedure for evaluating 
the distributions / in such a way of accounting for Pauli 
exclusion principle. This is easily done by imposing that 
a generic distribution / is given by the combined proba- 
bility of transmitting an electron with that of finding the 
final state empty. For instance the distribution of elas- 
tically transmitted electrons in the right-hand side lead 
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reads 



(f+ a ei) R ( E ') = F L (E')c(E')T aa (E',E') 



Ol'^OL 



E 

Ol'^O. 



U-if iE')} • (23) 



This means that the distribution in the right-hand 
side lead of electrons elastically transmitted at an en- 
ergy E' is proportional to the probability of transmit- 
ting an initially filled state in the left-hand side lead 
F L (E')T aa (E', E'), and to the probability that the final 
state is not already occupied [l — J2a'^u(f+ a in) R (E') — 

T, a >? a (f-in) R ( E ')] ■ This last probability is determined 
by two competing effects: inelastic transmission from the 
left-hand side lead {f+ c \ n ) R i and inelastic reflection from 

the right-hand side lead {f1 a in ) R . 

Similar equations can be derived for all the other dis- 
tributions in the right-hand side lead 

(f+in) R (E') - F L (E)c(E)T aa '(E,E>) x 



[1 



(f+ a el) R (E') ~ (f« a el ) R (E>) 



E u-i 

Ol'^O. 



) R (E)] , (24) 



(f- a ei) R (E') = Fji{E')d{E')R aa {E' , E') x 
[1 - E (f+in) R (E') - E (f-in) R (E')] , 



(25) 



Ol'j^OL 



ol' ^a. 



{n%) R {E') = F R (E)d(E)R aa (E, E') x 
[l-(fZi) R ( E ')-(f- a ei) R ( E ')~ 



{f-%) R {E') ~ E (f+%) R (E)] , (26) 



Ol' ^Ol 



where {f^ n ) R (E'), {f1 a el ) R \E') and (f?%) R (E') are 
the distributions for electron undergoing respectively to 
inelastic transmission, elastic reflection and inelastic re- 
flection, and we have introduced the reflection coefficient 
R a @ (E, E') = |r Q /3| 2 . Analogous expression for the left- 
hand side lead can be obtained by replacing R —> L, 



T a/3 rpapi = \f' a0 \ 2 and R al3 



R a P' = \ r ' a0 \ 2 in the 
equations (|23[1 through Q26J1. 

The normalization functions c(E) and d(E) are ob- 
tained by imposing charge conservation deep into the 
leads, which means that the sum of all the distributions 
/ must be equal to the equilibrium Fermi distribution of 
the particular lead 



F L (E) = (f™ l ) R (E) + (f™ l ) L (E) + 
J2(f+tn) R (E')+J2(f^n) L (E'), 



F R (E) = {f™ l )^(E) + {f™ l 

J2(f+in) L (E')+j2(f-in 



R {E)- 
R (E') 



(27) 



(28) 



The equations (|27|) and (|28|l together with those for the 
distributions (J2HJ , |2U , 123) and f^Hjl form a close set of 
equations, that can be solved for the / iteratively. These 
are evaluated numerically over a discrete energy mesh. 
Finally the current is obtained by a direct extension of 
the Landaucr formula, which takes into account the self- 
consistent electron distributions 



/ = 




(ft a el) R (E) +J2V%An) R (E) 



E (f+%l) L (E)+Y^(f+M) L W (29) 



III. RESULTS 

We are now in the position of investigating inelastic 
transport in atomic scaled systems. In particular we will 
consider three material systems. The first is a simple 
diatomic molecule. This will be used as a test case for 
explaining some general concepts of inelastic transport 
and for providing details about the convergence. Then 
we will investigate the case of the H2 molecule, for which 
experimental data are available [40| . Finally we will move 
our attention to the investigation of inelastic scattering 
in molecular spin-valves made from half-metallic leads. 
In all the calculations we consider temperatures consid- 
erably lower than the energy of the lowest of the phonon 
modes. This means that \k;0, 0, ...) is the only occupied 
incoming channel and that there cannot be a net phonon 
absorption. 



A. Fixing the phonon cutoff 

The first system that we consider is a diatomic 
molecule sandwiched between two identical mono- 
dimensional leads. These are described by a single orbital 
nearest neighbor tight-binding model with on-site energy 
«l = £r=-10 eV and hopping parameter /3l = /3r=2 eV 
(the total bandwidth is 8 eV). The Fermi level of the 
leads is taken at Ep=-10 eV, which corresponds to half- 
filling. The tight-binding parameters of the diatomic 
molecule are taken from reference [Tfjj and are cm=- 
10 eV, w°=0.15 eV, i 4j =0.2 eV, ay =6.1 eV/A. The 
masses of the atoms in the molecule are 13 a.m.u. (each 
site corresponds to a C-H group) and the spring constants 
are Ki cad -M=180 eV/A 2 , and K M = 90 eV/A 2 , where 
the first is the one between the leads and the molecule 
and the second is the intra-molecular one. 

These parameters lead to two longitudinal phonons 
(the only ones considered here) of energies fiwi=0.24 eV 
and Sw2=0.33 eV. A schematic diagram of the energy 
level alignment is presented in figure [SJ The main goal 
of this section is to discuss how the energy mesh and 
the number of phonon states included in the calculation 
affect the results. 
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FIG. 5: Schematic representation of the energy level lineup of 
the diatomic molecule discussed in this section. The solid box 
represents the density of states of the mono-atomic leads, the 
solid horizontal lines denote the position of the bonding and 
antibonding states of the free molecule and the thin horizontal 
line is the position of the Fermi level £"f- 



The general strategy for testing the convergence of 
the calculation is as follows. First one investigate how 
the transmission coefficient changes with the number of 
phonon states iVph included in the calculation. This is 
done with a reasonably fine energy mesh in order to have 
a good energy resolution. Then, once N p ^ is fixed, the en- 
ergy meshed is varied and the convergence of the current 
is controlled. Note that the energy mesh and the cutoff 
over the number of phonon states have two deeply differ- 
ent physical meanings. The energy mesh merely controls 
how accurate is the energy resolution, and typically a 
resolution of the order of 1 meV is sufficient. In contrast 
the number of phonon states essentially establishes the 
degree of accuracy in describing multiple phonon scatter- 
ing. 

In figure H3 we show the total transmission coefficient 
(elastic plus inelastic) at zero bias as a function of energy 
for different numbers of phonon states \rii, n-i) included in 
the calculation. It is clear that as N p ^ increases the trans- 
mission coefficient becomes more broadly distributed in 
energy and the resonances are largely smeared out. When 
no phonons are present T(E) has only two sharp reso- 
nances, located at the bonding and anti-bonding ener- 
gies of the molecule, with a broadening proportional to 
the electronic coupling with the leads. The inclusion of 
phonon absorption/emission opens new scattering chan- 
nel where phonon-assisted hopping to the molecular level 
becomes possible. This broadens the spectral range of 
T(E) and shifts the transmission peaks. In particular 
the resonance at the bonding state shifts almost rigidly 
of about 0.15 eV towards lower energies, while that at 
the antibonding level gets substantially smeared. 

This behavior is somehow expected since in our work- 
ing conditions of low temperature electrons cannot be 
transmitted with a net energy gain (no net phonon ab- 
sorption is allowed). Thus most of the spectral modi- 
fications due to electron-phonon scattering occur above 
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FIG. 6: Total transmission coefficient (elastic plus inelastic) as a 
function of the energy obtained for an increasing number of phonon 
states TVph. The convergence is obtained for 24 states (solid line), 
that is a maximum occupation of n a =5 for each phonon mode. 



the Fermi level. If we now focus our attention on the 
transmission coefficient for N p ^=24 (the higher occupa- 
tion allowed is n Q =5 for both the phonon modes) we can 
associate the first peak at about -10.25 eV with resonant 
elastic transport through the bonding state, the second 
just below -10.0 eV with phonon-assisted resonant trans- 
port through the bonding state via the first phonon mode 
(fkui—0.2A eV), the third at about -9.92 eV with phonon- 
assisted resonant transport through the bonding state via 
the second phonon mode (fia>2=0.33 eV), and so on. Note 
that peaks in the transmission coefficient appearing at 
rather high energies are the result of resonant transport 
through either the bonding or the antibonding level via 
multiple phonon emission. 

Note also that upon increasing the number of phonon 
states T(E) is only weakly affected for low energies, but 
changes drastically at high energies. This is because new 
scattering paths involving multiple phonon emission be- 
come available. However these high order multiple scat- 
tering events are increasingly improbable and one expects 
T(E) to saturate for JV p h large enough. As a measure of 
the convergence of T with the number of phonon modes 
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mission coefficients T m and T n calculated respectively for 
Nph^m and N p b=n 



M 



1 



M ^ 



T m {Ei) — T n (Ei)\ , 



(30) 



with the sum running over the energy mesh points. 
The values of A mn for the transmission coefficients 
of figure El are respectively A o , 3 =0.0242, A 3i8 =0.0121, 
A 845 =0.0065, A 15:24 =0.00 30 and A 2 4 : 3 5 =0.'0012. A 
value of A ~ 0.0015 is usually considered as a good level 
of convergence. 

The convergence of the transmission coefficients clearly 
guarantees the convergence of all the quantities related 
to them. Therefore after having determined iV p h we can 
then calculate the current using equation (|29|l . In figure 
|7| we show the elastic and inelastic contributions to the 
I-V characteristic and the total current (elastic plus in- 
elastic) as a function of A p h for different biases. A simple 
rule of thumbs for understanding the features of the I-V 
curve is that most of the current originates from the total 
transmission coefficient in an energy window comprised 
between Ep + V/2 and Ep — V/2. Note that this is only 
approximately valid when electron-phonon scattering is 
present. In fact the self-consistent electron distributions 
/ can provide some non-vanishing contributions to the 
current coming from transmission coefficients outside this 
energy window. 

First consider the case in which N p h=0, i.e. the purely 
elastic case. The I-V has only the elastic component 
and shows three different slopes, corresponding to three 
different resistances. For bias below 0.2 Volt the trans- 
mission is small, since no molecular state is comprised 
in the bias window. Then for voltages between 0.2 and 
0.4 Volt, both the bonding and the antibonding states 
of the molecule appear in the bias window, producing a 
large increase of the current. Finally for higher biases 
the current saturates since no new molecular states are 
available for transport. 

The main effects of switching on the electron-phonon 
interaction 7V p h 7^ is that the inelastic current starts 
to participate to the conductance. This is almost zero 
for biases up to about 0.2 Volt, which is roughly the 
energy of the first phonon mode, and then increase al- 
most linearly with bias. This opening of phonon-assisted 
transport channels appears in the I-V curve as a peak 
in the derivative of the differential conductance, a prop- 
erty which is used for detecting molecular vibrations from 
transport measurements [3^, [!(], EH ■ Note that phonon 
emission for biases smaller than the energy of the first 
phonon mode is certainly possible. However this process 
contributes little to the current, since the final state has 
high probability to be filled. Note also that the elas- 
tic current for small biases increases with respect to the 
purely elastic case. This is due to the contribution of 
multiple scattering, which effectively increases the value 
of the transmission coefficient around Ep (see figure EJl ■ 
The relative portions of elastic and inelastic current de- 
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FIG. 7: a) Elastic current for different phonon cutoffs, b) In- 
elastic current for different phonon cutoffs, c) Total current (elas- 
tic+inelastic) as a function of number of phonons. Each line rep- 
resents different voltages. At V= the total current is always 0. 



pend on the details of the Hamiltonian, however our anal- 
ysis demonstrate that a considerable number of phonon 
states must be included in the calculations for saturating 
the current. 

Figure [S] shows the same currents of figure this time 
obtained without imposing the Pauli exclusion principle 
through the distributions /. The most evident difference 
is that now the inelastic current does not vanish for V < 
0.2 Volt, i.e. for biases smaller than the energy of the 
first phonon mode. This happens since the occupation 
of the final state is not considered, thus transmission to 
final energies below E-p — V/2 (the chemical potential of 
the right-hand side lead) is possible. As a consequence 
the elastic current in the same bias region is reduced 
with respect to the previous case. Finally note that the 
currents calculated either imposing or not imposing the 
Pauli principle are similar for large biases, where the vast 
majority of the incoming states is filled, and the final 
states are empty. 



B. Transport through H2 molecules 

Platinum point contacts sandwiching H2 molecules 
represent the ideal system, where to investigate the ef- 
fects of electron-phonon interaction on the electron trans- 
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FIG. 8: a) Elastic current for different phonon cutoffs, calculated 
without considering the exclusion principle, b) Inelastic current 
for different phonon cutoffs, calculated without considering the 
exclusion principle, c) Total current (elastic+inelastic) calculated 
without considering the exclusion principle (elastic+inelastic) as a 
function of the number of phonons. Each line represents different 
voltages. At V= the total current is always 0. 



port. These contacts in fact are stable for long times, en- 
abling the measurement of low bias I-V characteristics 
under different stretching conditions [40|. The most im- 
portant aspect of the current /voltage curves is a rather 
sharp drop of the differential conductance Gr(V) = ^ 
from a value of approximately Go = 2e 2 /h for voltages 
close to 63 mVolt. This corresponds to a peak in the 
derivative of Gr, ^t 1 , and it is attributed to the onset 
of electron-phonon scatte ring to the lowest phonon mode 
available in the junction |40| . 

Despite the initial controversy recent oh initio 
calculations |43l 0] have convincingly proved that the 
transport in this system occurs through the antibonding 
state of the H2 molecule in the so called "bridge" con- 
figuration, i.e. with the molecular bond lying parallel to 
the transport direction. The antibonding state provides 
a high transmission channel and the differential conduc- 
tance is about 0.9 Go- Moreover ab initio calculations 
have also suggested that a transversal vibrational mode 
of the H2 molecule is responsible for the differential con- 
ductance drop at ~60 meV 0]. This last identification, 
although supported by the blue-shift of the differential 
conductance drop upon stretching, however is more dif- 



ficult to establish since the local arrangement of the Pt 
atoms at the tip is unknown (it is likely that the H2 
molecule bridges the Pt contacts in a zig-zag geometry 

In an attempt of understanding the main mechanism 
responsible for the drop of Gr(V) and to provide a fur- 
ther test to our scheme, we model a Pt-H2-Pt point 
contact as two semi- infinite monoatomic current /voltage 
leads sandwiching a H2 molecule. The leads are at half- 
filling with on-site energy eo=0 eV and hopping param- 
eter 7=5.0 eV. This guarantees a rather large (20 eV) 
bandwidth. The tight-binding parameters for the H2 
molecule are chosen in such a way of aligning the an- 
tibonding state at the Fermi level and to give a conduc- 
tance close to 0.97 Go for zero bias [ipj. We consider 
as relevant phonon modes only longitudinal modes and 
we tune the electron-phonon coupling in such a way of 
reproducing the experimental drop of the differential con- 
ductance. Note that at this level of approximation the 
choice of using either transverse or longitudinal modes is 
somehow immaterial. 
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FIG. 9: Differential conductance for the H2 molecule sandwiched 
between Pt leads. The solid line is for the full calculation where 
both the Pauli principle and electron-phonon interaction are in- 
cluded, the line labeled with "phonons" is for a the calculation done 
without considering the Pauli principle and the line labeled with 
"no phonons" is for purely elastic transport. In the inset the drop 
of the differential conductance as a function of the electron-phonon 
coupling a. The relative drop A is defined as the ratio between the 
differential conductance at V just above the first phonon energy 
(V=63 mV) and that at zero bias A = [Gr(V) - Gr(0)]/Gr(0). 

Figure |2| summarizes the main result of our calcula- 
tions. Clearly, when both electron-phonon interaction 
and Pauli principle are considered the differential con- 
ductance has a drop for voltages around 60 mV, which 
corresponds to the energy of the first longitudinal phonon 
mode at 63 meV (oscillation of the center of mass of 
H2). This drop of Gr is due to the enhancement of in- 
elastic back-scattering, which reduces the transmission 
of the elastic channel resonant at the antibonding state 
and it has the net effect of suppressing the current. In 
fact a careful analysis of the electron distributions / has 
revealed a peak in /£ in for energies corresponding to 
the emission of the first longitudinal phonon. This emis- 
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sion, leading to the enhanced inelastic back-scattering 
is Pauli suppressed for biases V such that eV < hio ll 
but it becomes possible as soon as the bias window is 
large enough to allow the phonon emission. Furthermore, 
within our simplified model we do not expect any addi- 
tional differential conductance drop for small bias, since 
the next longitudinal phonon mode has a rather high en- 
ergy (^430 meV). 

Note that such a drop of Gr for eV ~ Huii can- 
not be described unless both electron-phonon coupling 
and Pauli principle are considered. On the one hand, 
when the Pauli principle is neglected the inelastic back- 
scattering is always present and contribute to the current, 
with the net result that no conductance drop is found. 
In this case the conductance has a value similar to that 
obtained for V > 63 mV. On the other hand, if electron- 
phonon interaction is set to zero, no drop can be found 
since inelastic back-scattering is not included. The con- 
ductance now is similar to that calculated for V < 63 mV 
and almost independent from the bias. 

Importantly the magnitude of the differential conduc- 
tance drop depends on the strength of the electron- 
phonon coupling, and it grows as the coupling is en- 
hanced. Such a behavior is demonstrated in the inset 
of figure where [Gr(V) — Gr(0)]/Gr(0) is plotted as a 
function of the coupling a. Interestingly we found that 
A a a 2 . Since the relevant longitudinal phonon mode is 
that of the center of mass of H2 oscillating between the 
two Pt contacts, the only a determining the transport is 
that of the Pt-H bond. This is mainly an sscr bond |43[ . 
which roughly scales as 1/cP with the separation d be- 
tween Pt and H Such a scaling leads to A oc l/rf 6 , 
and therefore we predict for longitudinal modes a very 
drastic dependence of the differential conductance drop 
as a function of the point contacts separation. 



C. Molecular Spin-valves 

Finally we investigate the effects of inelastic electron- 
phonon scattering on the transport properties and the 
GMR of a model spin-valve formed by half-metal cur- 
rent/voltage probes sandwiching an organic molecule. In 
particular we consider a Mp half- metal [27|. in which 
the Fermi level cuts through the minority band (see fig- 
ure In our model calculation this is obtained by a 
linear chain of hydrogenic atoms whose on-site energies 
are spin-split. In particular we consider = — 14 eV for 
the majority (|) and = —6 eV for the minority (|) spin 
bands. For both the spin-bands the hopping parameters 
are /3l = /3r = — 2 cV, which give a bandwidth and an 
exchange of 8 eV. With this choice of parameters the ma- 
jority upper band-edge coincide with the lower minority 
band-edge at E=-10 eV. 

The molecules are linear chains of various lengths (from 
one to three atoms) made by atoms of mass 13 a.m.u. 
(C-H groups) and characterized by the hopping inte- 
gral tij — 0.2 eV. The on-site energy of the molecule 



is always close to the position of the Fermi level giv- 
ing resonant transport through the molecular states. We 
have then investigated two possible situations. In the 
first one (see figure HUk) £m = —8.2 eV and the Fermi 
level Ep — —8.0 eV is well within the minority band. 
In the second one Ep = em = —9.7 eV is close to the 
lower minority band edge ffigure lTUb ). The hopping pa- 
rameters between the molecule and the leads and the 
electron-phonon couplings are the same as reference |16| , 
i.e. w°=0.15 eV, a = 6.1 eV/A, K lcad ^ M = 180 eV/A 2 
and Km = 90 eV/A 2 . This leads to phonon modes with 
the following energies: i) one atom chain, fiwi=0.34 eV; 
ii) two-atom chain ftu>i=0.24 eV, hoj2=0.33 eV; iii) three- 
atom chain Huji=0.17 eV, ?iu;2=0.30 eV and So;3=0.34 eV. 




FIG. 10: Schematic representation of density of states for the 
model Mpi half-metal and the two configurations studied: a) Fermi 
level well within the minority band, b) Fermi level close to the 
minority lower band-edge. 

In investigating the spin-transport we consider the two 
spin fluids approximation where spin-flip scattering is ne- 
glected and the two spin sub-bands conduct in paral- 
lel. This is a good approximation for C-based organic 
molecules since both spin-orbit and hyperfine interac- 
tions are weak. The GMR is then defined as usual as 

MR = , (31) 

Jp + Jap 

where Ip is the current of the spin- valve when the mag- 
netization vectors of the two leads are parallel to each 
other and Jap that for the antiparallel alignment. Within 
the two spin fluid model these read Ip — P' + P* and 
Zap = I + 1 , with I a " the spin current and the index 
a ((3) indicating the spin direction in the left-hand side 
(right-hand side) lead. Thus for instance P* is the spin 
current for the antiparallel configuration where electrons 
propagate in the majority band in the left-hand side lead 
and in the minority in the right-hand side lead. 
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Let us consider first the case when em = —8.2 eV and 
Ep = —8.0 eV. In this situation the current is not neg- 
ligible only for the minority spins in the parallel config- 
uration 1^ at least for biases lower than 3.6 Volt, after 
which also P^ gives contributions. In fact the particular 
density of states of our model half-metal and the rela- 
tively large Fermi level, result in the fact that for all the 
other configurations the density of states (DOS) of one of 
the two leads vanishes within the relevant bias window. 
The current is therefore completely suppressed by the 
lack of available states in the leads, and the GMR is al- 
ways 1. Note that in principle multiple phonon emission 
might result in finite current also for other configurations, 
however this is always entirely suppressed by the Pauli 
exclusion principle. 
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FIG. 11: I-V characteristic for the model spin-valve with Ep = 
€M = —8.2 eV described in the text. The molecules are formed re- 
spectively from one atom (upper panel), two atoms (middle panel) 
and three atoms (lower panel). The only non-negligable contribu- 
tion to the current is from the minority electrons in the parallel 
configuration 1^*. 

The calculated I-V for the parallel case (up to 1 Volt) 
is presented in figure HTI for the three atomic chains inves- 
tigated. The most notable feature is the drastic suppres- 
sion of the elastic current as the chain length increases, 
and the fact that the inelastic component is only weakly 
affected. This is somehow expected since as the molecules 
becomes longer the direct tunneling probability decays 
and the current is dominated by the inelastic component. 
Interestingly, for this choice of parameters, reproducing 
(CH)„ molecules, the inelastic current dominates already 
for rather short lengths (n=3). 

The second case, when Ep — e M = — 9.7 eV is some- 
how different, since the Fermi level cuts much closer to 



the minority lower band-edge. In this situation the half- 
metal DOS suppresses the current only at low bias and 
one expects transport in the antiparallel configuration, 
due to the 1^ channel for biases as small as 0.6 eV. In 
fact at such biases finite DOS appears in the bias window 
for both the leads. This does not happen for the P* and 
P' contributions, which vanish at all biases. Therefore 
one has Ip — 1*^ and Jap = P^ and the GMR is in 
general dependent on the bias. 




V(W) 

FIG. 12: I-V curves for the parallel and antiparallel configura- 
tions of a half-metal molecular spin- valve made from atomic chains 
containing (a) one, (b) two and (c) three atoms. In this case 
Ep = em = —9.7 eV, Ip = 1^ and I\p = I'*. In the insets 
the corresponding GMR as a function of bias. 

Figure ^| shows the I-V curves for both the paral- 
lel and antiparallel configurations of the spin-valve and 
the corresponding GMR, for molecules made respectively 
from one, two and three atoms. Let us consider first 
Ip = p±. The elastic component of the current dom- 
inates for voltages up to approximately 1.5 Volt, after 
which it gets completely suppressed. This suppression 
originates from the fact that the "effective" bias window 
(where both the leads have finite DOS) remains constant 
for V > 0.6 Volt, but it moves increasingly far from 
the energy levels of the molecular states. This drasti- 
cally reduces the elastic transmission coefficient with a 
strong decrease of the current. The same effects some- 
how occurs for the inelastic component of the current al- 
though phonon emission allows transmission over a much 
wider energy window. This makes the inelastic compo- 
nent dominating for V > 1 Volt, regardless of the chain 
length. 

In contrast in the antiparallel configuration no current 
can flow until the bias is large enough that the DOS of 
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the two leads overlap for some energies. This happens for 
V ~0.5 Volt for the elastic component and V ~1.0 Volt 
for the inelastic one. In this case we do not expect a 
decrease in the current for increasing bias, since the bias 
window is not cut by any band edge. As in the previous 
case the inelastic component becomes progressively more 
dominant as the chains gets longer and it is expected to 
dominate completely for very long chains. The result- 
ing GMR is strongly bias dependent. In particular, for 
the present choice of parameters, it changes sign from 
positive at small bias to negative a high bias. 
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FIG. 13: Elastic I-V curves for the parallel and antiparallel con- 
figurations of a half-metal molecular spin-valve made from atomic 
chains containing (a) one, (b) two and (c) three atoms. In these 
calculations no electron-phonon interaction is considered. Ep = 
em = —9.7 eV and Ip = i"H and I\p = JTi, l n the insets the 
corresponding GMR as a function of bias. 

As a comparison in figure (|13J) we present results for 
the same systems obtained without electron-phonon in- 
teraction. In this case there is no sustantial reduction of 
the elastic current as the molecular chain becomes longer, 
and the transport is always resonant through the molecu- 
lar levels. However in absence of the inelastic component, 
the total current for the parallel alignment sharply decays 
for biases in excess of 1 Volt. As a consequence the GMR 
sharply changes from 1 to -1 as the bias increases. 

Generally, inelastic scattering increases the spectral 
range (the energy region where the total transmission co- 



efficient does not vanish) contributing to the current (see 
figure El- In this case of half- metallic leads, the effect 
is to change the GMR as a function of bias. However if 
the leads are made from transition metals, then the effect 
of electron-phonon scattering is expected to be a change 
in the spin-polarization of the current pfij ] . This is con- 
nected to the increase of the spectral range of the trans- 
mission coefficient. A large spectral range means that 
a large portion of the ferromagnetic DOS contributes to 
the current. However in an ordinary transition metal the 
spin-polarization of the DOS is large at around Ep but 
decreases for both higher and lower energies. At higher 
energies this is due to the fact that there is little contribu- 
tion from the d electrons for both spin directions, while 
at lower energy a large g?-DOS is present for both spin 
directions. We therefore expect that electron-phonon in- 
teraction will produce a general suppression of the GMR. 



IV. CONCLUSIONS 

We have demonstrated the role of inelastic electron- 
phonon scattering on the transport properties of differ- 
ent nanoscale devices. This is treated using the multi- 
channel scheme introduced by Bonfia and Trugman with 
the explicit inclusion of the Pauli exclusion principle. 
Such method allows us to map a scattering problem with 
inelastic effects, onto a fictitious problem where only elas- 
tic scattering events are present. We have then given a 
full description of the method, within our surface Green's 
function formalism, and described the criterion of conver- 
gence over the number and nature of phonon modes. 

Two specific material systems have been investigated. 
The first are Pt point contacts sandwiching a H2 
molecule. In this case the differential conductance shows 
a sudden drop for a bias corresponding to the energy of 
the first H2 longitudinal phonon mode, which we have 
explained as an enhancement of inelastic backscattering. 
Note that such a drop cannot be found unless the Pauli 
principle is considered. 

Then we have investigated two types of spin- valves, 
obtained by contacting with half-metals a linear chain of 
CH groups. In this case the GMR depends on the specific 
alignment of the Fermi level with respect to the minority 
lower band-edge, and in general it is strongly bias depen- 
dent. The presence of electron-phonon scattering largely 
increases the current in the parallel state at high voltages 
and makes the change in GMR upon bias smoother with 
respect to the case where electron-phonons interaction is 
not considered. 
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